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Abstract 

We investigate effects of strong interactions on tlie electromagnetic dissociation 
of nuclei in heavy ion collisions. 

We start from the eikonal approach to the equivalent photon method to describe 
the electromagnetic contributions to the cross section of peripheral collisions. We 
summarize some results of this approach and characterize recent experiments in 
"universal plots". 

In the second part of this work, we give a straight forward method to include 
transitions induced by the strong interactions between the ions. We introduce dif- 
ferent methods to obtain the nuclear transition potential, and study the behavior 
of the resulting nuclear contributions. 



1 Introduction 



The Coulomb field of a fast moving heavy ion is an intense source of quasi real photons 
ll], 0, 1^. With raising projectile energy this photon spectrum becomes harder and there 
exist a variety of applications in nuclear- as well as astrophysics, like the excitation of 
giant resonances or the Coulomb dissociation of the ions [Q. 

The advantage of using heavy ions to study this processes are the large cross section, 
which make the Coulomb dissociation an interesting tool to study even multi-phonon 
resonances ||^, |^ or the astrophysical S-factor for reactions like Be{p,'-f)^B [Q. In the 
latter case phase-space considerations and the large number of equivalent photons highly 
favor the Coulomb dissociation compared to the direct measurement of the radiative 
capture reaction |[. 

A powerful but simple tool to describe the Coulomb dissociation of heavy ions is the 
equivalent photon method . On the other hand this simple approach does not take into 
account excitations due to the strong interaction. Therefore it is interesting and necessary 
to study these effects. Using the hydro-dynamical model of Bohr and Tassie for the nuclear 
transition densities or potentials P], one can connect the nuclear interaction to the optical 
potential and the deformation parameter and therefore to the electromagnetic transition 
matrixelements. This approach allows us to study these contributions in various cases of 
physical interest. 

2 Equivalent photon spectrum 
in a Glauber approach 

In the following we consider a situation where a projectile (particle 1) with charge Zi 
excites a target (particle 2) with charge Z2 from a state |/j = 0, Mi = 0) to the excited 
state I J, M) in a peripheral collision via the exchange of one quasi real photon of energy 
= huj. The velocity /3 = ^ of the projectile is in z direction and we use 7 = (1 — ^ . 
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The scattering angle of the projectile is denoted by 6. The inelastic scattering amplitude 
in the eikonal approximation f{6) is than given by |jlO|, |lT|: 



f{e) = ^Jd'r d\' e'^-^^ e'^^'\^f{T')\VUr.r')mr')). (1) 

where q denotes the momentum transfer and k is the wavenumber of the incoming pro- 
jectile. We define r to be the separation of the centers of mass of the two nuclei, and r ' 
to be the intrinsic coordinate of the target nucleus. 

The Glauber phase x(^) is the sum of a Coulomb and a nuclear part. If we assume 
the 'sharp cut-off' model for the nuclear part 

1 1 iih> R 

where i?j = 1.2 fm A^/^, (z = 1, 2) is the radius of the i'th nucleus, the Coulomb phase is 
given by 

Xcih) = 2 r, ln{kh). (3) 

Here rj = ^^^^^^ is the Sommerfeld parameter. 

Using standard methods 0] and defining the reduced electromagnetic transition prob- 
abilities B{ttI) = B{ttI, li If) ||12| we can write the cross section in the form 



dn-\J ^""^^^^^^In- 



In the above formula a^\E^) is the photo-nuclear absorption cross section for a given 
multipolarity 



where priyj) is the final state density. 

Under these assumptions the equivalent photon number per unit solid angle for a given 



multipolarity is given by: 



M = 72 /a;fcy /[(2/ + l)!!P ^ 
dn '"UV (2vr)3(/ + l)^ 



G 



1^ ' 



.Irr. 1 ^ 



fi^(gT)r, (6) 
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Here the transverse momentum transfer ^ kO is the only quantity, which depends on 
the scattering angle. Defining the variables 

K 7f 

'^'^'^^ ^ ~kR ' ^^''^ '^^^'iiff^ C^) 

can be written as: 



J I Odiff 



In the case of the classic trajectory of the projectile {rj oo), one can use the saddle 
points method to solve the integral in eq. (|8D[pT|: 

= R'^ KU^y) , (9) 

where the variables y = ^ and = ^7r{m + 1) + 2ri {in [jf^ — l) where introduced. 

Expression (eq. ^ for Vtrn mainly depends on two parameters, the interaction strength 
7] and the adiabaticity parameter ^, which is the ratio of impact time to transition time 
in a grazing collision. 

It is interesting to compare |2r7^p as a function of j- for fixed ^ but different i] as 
shown in the 'universal plots' in figs. (1) and (2). The thick solid line shows the semi 
classical limit according to eq. (^) while the dashed and dotted lines are the solutions of 
eq d^) for different values of 77. At ^ = 1 the semi classical limit drops to zero while the 
full solutions extend to higher angles. At small angles the full solution shows diffraction 
effects, but for a large interval of angles the full solution for sufficiently large r/ is a fast 
oscillating function around the semi classical limit (dotted line). Only for small rj (dashed 
line) the full solution differs substantially from the semi classical limit. 

For the above reason the semi classical solution for the inelastic scattering amplitude 
(eq. P provides in many relevant cases already a good description of the experimental sit- 
uation. Fig. (3) classifies some existing experiments in a ^-?7-plot where the experiments 
well described by the semi classical solution are at large r] and small ^. 



To complete the discussion of the sharp cut-off model we want to give the expression 
for the total cross section, which is obtained by integrating eq. (|^) over the solid angle 
and using the closure relation for the Bessel functions. One finds: 

atot = 27: [^-^j ^^j ^(^^j \G^im{^)\'B{7rl)R' x dx KU^x). (10) 

Notice, that this expression is independent of r]. 



3 The nuclear excitation part 

The sharp cut-off model, as it was introduced in the last section, takes nuclear effects into 
account only in a very simplified way. A more realistic model for the nuclear contributions 
will influences the inelastic scattering amplitude (eq. |I|) in two different ways. First it 
leads to an additional transition potential in eq. ([^) of the form {^f{f')\Vj\f\^i{f')), where 
Vn is the nuclear potential between the two nuclei. In addition it will change the nuclear 
phase due to 

XNib) = -j-J^ V^{yWT^)dz, (11) 

which, in contrast to the sharp cut-off model, now has a finite real and imaginary part. 
This leads to a "smoother" cut-off than in the sharp cut-off model and a modification 
of the Coulomb potential Vc due to the penetration of the two charge distributions. In 
the following we will neglect the effect on the Coulomb potential but keep eq. (|n|), and 
concentrate on the nuclear induced transitions. 

There exist different approaches to obtain the nuclear transition potential (eq. ^ 
($j(r ')|VAr|<l'j(r ')). In this work we assume a collective nuclear model. 

First we want to describe the approach via the transition densities. In this case the 



nuclear transition potential can be calculated using the folding formalism [[13| : 

($/(r-")|lV(r,r-")|$.(r-")) = (t^7v(^)) Pi{r-f') {I , M\p2{f')\0,0) . (12) 
Here pi, {i = 1,2) denotes the charge density of the i'th nucleus. For the energy de- 



pendent tATAT-matrix see ref. [13|. To determine the transition density matrix elements 



{I , M\p2{f')\0,0) we assume the hydro-dynamical model of Bohr and Tassie [Q: 

(/,M|p,(f')|0,0) = ^gti.'C-"^ y;„(f 0. (13) 

The deformation parameter j3i is related to the electromagnetic transition matrixelements 
M{El,m) by: 

M{El,m) = VWTiml'' / dr r^^alr) 6M,m (14) 

Jo 

Under these assumptions the inelastic scattering amplitude (eq. |1|) can be written as: 

„-m+l I. 

fi0)= E^im ^— e-^'"" (/,M|M(7r/,m)|0,0) 

/o°° 6ci6e*(^^(^)+'^'^(^»J^(gr&) {r^''"(6) + r^"'(6)} , (15) 

where a is the azimuthal scattering angle. Since we neglect the corrections of the Coulomb 
part due to the penetration of the two particles, T'^"^ is still assumed to be given by: 



7 c \f3J 



(16) 



For the sake of mathematical simplicity we consider a Gaussian shaped charge distri- 
bution for the two nuclei [l^ to calculate F]!^'": 



p,(r) = p,(0)e ^ ; p,(0) = O.OSS/m-^e^^ . (17) 

We have used the notation R^. = 2aiRi, where Oj is the diffuseness parameter of the 

i'th nucleus. This choice for the charge distribution allows us to perform most of the 
calculations in an analytical way. 

Under the above assumptions the nuclear phase (eq. |Tl]) is given by 

XN{b) = (tME)) Pi(0) P2(0) 1%?^} e-''/''o^ (18) 

where Rq = Rq_^ + Rq^. To take into account the singularity Xc{b = 0) (eq. |^), we modify 



the Coulomb phase |11 



i 



[ (4 



Xc(b) = 2r, { In(kb) + | . (19) 



In the above formula Ei denotes the exponential integral. 

Calculating the nuclear transition density as given in eq. (pTSD and using eq. ([T^) we 
find for r]^'"(6) : 

Tf'^ih) = 0. (20) 

The Gauss parameterization for the charge densities is a relatively good approximation 
for light nuclei, but fails for heavy nuclei. Starting from more realistic densities immedi- 
ately makes it necessary to solve the folding integral (eq. |12D numerically. To avoid these 
difficulties one can start from optical potential rather then their corresponding densities. 
This formalism is described in ref. |]T3|. 

For deformed optical potentials of the form VN{R{0,(f)),r), where R(9,(l)) = Ro{l + 
Y^Xfi^^Xfi^Xfii^j (p)): that depend only on the distance {RiO^cj)) — r), the nuclear transition 



potential (/, M| V^|0, 0) is of the form p 



(I,M\V,M = (21) 



In both, the Tassie model and the Bohr Mottelson model [|T6|, 13 1 is than given by 



3 Z2Ri V e2 
Using eqn. (|21|, ^ |15|, and § we find for rf'™(6) : 



with = 6^ + 2;^. 

This result can be compared to the Distorted Wave Born Approximation (DWBA), 
see for example ref. |T^. In this formalism the nuclear part of the inelastic scattering 
amplitude for a transition from a state |/j = 0, Mj = 0) to the exited state |J, M) is given 
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by 

^'■^'^^ ^ ^ hv Y V2/+ 1 A;A;' 

^^^(A' /x' / M|A 0)(A' / 0|A 0) 1av(^>0), (24) 

where A; {k') is the wave number of the incoming (outgoing) wave. The radial matrixele- 
ment Ri'x has the form |]T9 



Ri, = /#.A,(t.V) ''^"'y-'-' />(^r). (25) 

For small scattering angles ^ and large values for the angular momentum of the partial 
waves A, A' one can show that eq. (p^) is equivalent to the derived Glauber approach 



(eqn. |T|, ||) 



This connection of the Glauber model to the DWBA can be seen as following. In eq. 
5]) the radial Coulomb wave function fx in the WKB approximation is given by [|l^, |20 

fx^ e'^^{k^/g{r))^sin(l> (26) 



where 



and tq is defined by gijo) = 0. 

Assuming R{,,^ ^ i?^^ with A = |(A + A') for sufficiently large values of A one can 



follow the calculations from ref. |20|, neglecting the rapidly scillating terms, and obtain 



4 J — OO UjRq 



where we have used e^^^^ = e*^^^''-* and neglected the Coulomb part in eq. (pTf), which 
leads to straight line trajectories instead of Rutherford trajectories. The A-dependence of 
i?^^ is contained in the impact parameter dependence of eq. (pS]), since one has A + | ~ kb. 
Making use of the approximation 



/ 2A' -I- 1 

^'-^(^,0) ^ J^^JMikbe), (29) 



in the limit of ^ 1 and replace J2x ^ k J db on has 

-ik 2 



' ' nv yvr J 

E^^-^'(A' - M/M|AO)(A'0/0|AO)i r W ^ 

4 J —oo 



dRQ 

Here we have also identified e**^'^^"'"'^^'-' ~ e^*""^ = e*^^^''-*. Finally we make use of the 



approximation ||19| 



X:^"-"'(A' - MIM\X 0)(A' / 0|A 0) ^ r^y_il_y,^,(|, 0), (31) 

for (/ + M) even and zero otherwise, and obtain for optical potentials Vn{R{0, 0), r) that 
depend only on {R{9, 0) — r) 



fN,i,M (^)-(-^) f,^ V2/ + li '^JM{khe)j_Jze^ Y, 

(32) 

This result is identical to eqn. (|15|, ^3]) for a surface peaked nuclear interaction {YiMidf, 0) - 
^/M(f )0)) apart from an overall phase, which doesn't infiuence the cross section. 

We want to stress, that we made essential use of the assumption, that the angular 
momentum of the contributing partial waves A and A' are large, and that the cross section 
is limited to small scattering angles, which both are sufficiently fulfilled for most heavy 
ion collisions. 

To apply eqn. (^, ^3]) for the nuclear contribution we want to assume a Woods-Saxon 



IM\ 



IT 



parameterization for the nuclear potential Vn = U{r) + iW{r) where 



lit 



U{r) = -Uo i^e +lj ; W{r) = -Wq i^e +lj . (33) 

In fig (4) we show the results with the above formalism for the Coulomb dissociation 
of ^B. This reaction has recently been studied experimentally in order to determine 



the astrophysical S-factor of the '^B{p,'~f)^B reaction |T^, The figure shows the cross 
section for a 51.9MeV/nucleon ^B beam on a "^^^Ph target. The excitation energy is 
assumed to be 1.2MeV. The solid line shows the total cross section {E1+E2), the dashed 
line the sum of the Coulomb and nuclear contribution for / = 2 {E2c+n), and the dotted 



line gives the E2n contribution only. We used Si-j{El) = 20eVb and Sij{E2) = lOmeVb. 
For the Woods-Saxon potential we assumed the following parameters |P2|| , 

Uo = 50MeV ; Wo = 57MeV ; 
Ru = Rw = 8.5/m ; au = aw = 0.8/m. (34) 

One should bare in mind that the use of the collective model for this —^^ Be + 
p continuum transition is not well justified. A diffraction dissociation approach would 
certainly be more appropriate, see e.g. ref. |^ . 

A more schematic model for the optical potential V/v = t^Ar + iW^ is given by the 
"ramp" potential, 

' -f/o \ir <R'^ 

f^(^) = -t^^^v -r) \iR-v<r<R^ . (35) 
if r > i?+ 

where R^ = Ru±^. The imaginary part Wj^ is given in a similar way. 

Fig. (5) compares the nuclear phase |e*^^*^'')p and |r^(6)p for Woods-Saxon (solid 
line) and the "ramp" potential (dashed line). The functions where calculated at = 2.0°, 
for Z = 2, m = 0, and are plotted as a function of hjR. For the Woods-Saxon potential 
we used = 50MeV, Wo = 57MeV, Ru = Rw = 8.5/m, and au = aw = 0.8 and the 
same parameters for the "ramp" potential but Au = A^y = 4/m. One finds, that for an 
adequate choice of parameters both models give similar results. 

Because the product \e'^^^''^T^^{b)\ is a good measure for the nuclear transition strength, 
we want to deduce some basic features of this quantity. The magnitude of the real part 



Uo is directly proportional to the magnitude of |r}^(6)|, while the imaginary part W, 



N 



mainly infiuences the shape of |e*'^^('')p. Namely for decreasing Wo the contributions for 
b < R — Aw in eq. (|T5|) becomes important. This leads to an increasing nuclear contri- 
bution, while the oscillations of the Coulomb contribution around the semi classical limit 
(figs. 1 and 2) decrease. For the radii a choice like Ru > Rw "will lead to a larger overlap 
in fig. (5) and therefore to a larger value for |e*^^*-^''r^(6)|. 



9 



This shows, that the magnitude of the nuclear contribution can vary over a wide range 
if one uses different sets of parameters for the optical potential. 

For sufficiently large Wq and small A one can deduce a very simple expression for the 
"ramp" model. For simplicity we assume Ru = Rw = R and Au = Aw = A. The 
nuclear phase is then approximately given by 

'0 iibKR^j 
'-f- iiRu<h<Ri , (36) 
1 iib> R-^ 

Using this approximations, we find for e^^^'^^'Tfl'^ih): 

Yi^{-,0)—sm I dj^i, (37) 

2 u \ V J 



for R^ <b< R]} and zero otherwise. In eq. (|37D we used F;m(^f,0) ~ ^jn(f ,0) and we 



have defined z^ax = y R^^ — b'^- 

Expanding the sinus for small arguments and introducing the variables D = {Ri 
R2)/R, X = b/R, and 6 = A/R one can write eq. 



(x - 1 + 5)^(1 + ^)2-x2 6j,i. (38) 
To see the influence of the parameters S and D Fig. (6) shows |/jv^^(^)P obtained 



from eq. (l37|) as a function of 9/Ogr for different values of 6 (solid line 6 = 0.5, dashed line 
5 = 0.1 and dotted line 6 = 0.02). We used D = 0.5, ^ = 0.1 and t] = 9. Fig (7) shows 
the same for fixed 6 = 0.5 but varying D (solid line D = 0.5, dashed line D = 0.3 and 
dotted line D = 0). 

We want to stress, that r^''"(6) is nearly independent from the excitation energy lj, 
while r^'''"(6) decreases with increasing excitation energy due to the behavior of Km 
(eq. HD. 



10 



Finally we want to discuss the limit of the sharp cut-off model for the nucleus in the 
adiabatic limit {k' = k). In this limit eq. (pSf ) becomes 



RUk' = k)= e^^-^'>- dz (39) 
and one can derive a relation between this radial matrixelement and the nuclear phase 

- - X V- '^"^ 

where e^^^''^ is like the optical potential Vn{R{0^ 4>)^f) ^ function of the distance {R{9, (f)) — 
r). Using eq. (j|) for the nuclear phase, one finds for the sharp cut-off model 

/7v,/,m(^) = e-''''^ (3ikR2R e^^^(«) FiM(^, 0) JMigrR)- (41) 
Using the asymptotic expansion for the Bessel function for large qtR one finds the typical 



oscillatory behavior for the nuclear contribution as predicted by Blair [1^ . 

Fig. 8 compares the results for the E2 nuclear cross section for the different models 
for the same experiment as in Fig. 4. {51. 9 MeV/nucleon beam on a "^^^Pb target). 
The solid line shows the result according to eq. ( |23D for a Woods-Saxon optical potential 
with the parameters as given in (p^. The dashed line shows the result for the ramp 
model with Uo = 50MeV, Wq = 57MeV; R^ = Rw = 8.5fm and Au = Aw = 4:fm 
and the dotted line shows the result for the sharp cut-off model (eq. ^l]). The sharp 
cut-off result shows the typical Frauenhofer diffraction while both other models show a 
smoother nuclear contribution and a phase shift at higher scattering angles as compared 
to the sharp cut-off model due to the diffuse edge of the optical potential. 



4 Conclusions 

In the first part we discussed the pure Coulomb contribution in heavy ion collisions using 
the Glauber approach. Using so called "universal plots", it is possible, to characterize 
this Coulomb contributions to peripheral heavy ion collisions for different experiments. 
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The relevant parameters are the interaction strength r), the adiabaticity parameter ^, and 
9/9gr- It turns out, that for many experiments a semi-classical approach to the equivalent 
photon method is already a sufficiently good description of the data. 

In the second part, we could include the nuclear excitation in a straight forward way. 
We used the folding formalism and the optical potential model, to determine the nuclear 
transition potential. The latter result was compared to the DWBA formalism. The 
introduction of the simplified "ramp" model makes it possible to clarify the role of the 
parameters of the used optical potentials. 

A comparison of the sharp cut-off model results with results obtained from different 
parameterisations of the optical potential shows, that the nuclear contribution can be 
described within the discussed frameworks. Therefore the introduced model for the nu- 
clear contribution makes it possible to estimate the nuclear effects on the electromagnetic 
Excitation and Dissociation of nuclei. 
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Figure Captions 

Figure 1: |2?7^|^ for = 0.5, r] = 20 (dashed) and r] = 2 (dotted) upper picture and for 
^ = 0.1, ?7 = 3 (dashed) and t] = 1.2 (dotted) lower one. The solid line shows the 
result for the semi classical expression (eq. 

Figure 2: Same as fig. (1) but for m = 2 

Figure 3: Some experiments characterized in a ^-r^-plot. 

Figure 4: Cross section for the Coulomb dissociation of 51.9MeV/nucleon beam on a 
^°^P6 target. The solid line shows the total cross section, while the dashed gives the 
E2c+N, and the lower line only the E2]\f contribution 

Figure 5: |r^^°(6)|^ (left hand side) and |e*^'^*^^)|^ (right hand side) as a function of The 
solid and the dashed line belong to the Woods-Saxon model and the simplified ramp 
model respectively. 

Fi gure 6: |/iv'^^(^)P as a function of 6/6gr for different values of S (solid line S = 0.5, dashed 
line 6 = 0.1 and dotted line 5 = 0.02) 
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Fi function of 6/6gr for different values of D (solid line D = 0.5, dashed 

line D = 0.3 and dotted line D = 0). 

Figure 8: E2 nuclear cross section for the various models with the same experiment as in fig. 4. 

The solid line shows the result according to eq. (|2^) for a Woods-Saxon Potential, 
the dashed line for the ramp model and the dotted line for the sharp cut-off model 
(eq. ^l]). The parameters used for each potential are given in the text. 
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